Dynamics of scalar perturbations in f{R,T) gravity 
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In the context of f{R, T) theories of gravity, we study the evolution of scalar cosmological pertur- 
bations in the metric formalism. According to restrictions on the background evolution, a specific 
model within these theories is assumed in order to guarantee the standard continuity equation. Us- 
ing a completely general procedure, we find the complete set of differential equations for the matter 
, density perturbations. In the case of sub-Hubble modes, the expression reduces to a second-order 

D ■ equation which is compared with the standard (quasi-static) equation used in the literature. We 

Mh ' show that for general f{R, T) Lagrangians the quasi-static approximation yields to very different 

\ results from the ones derived in the frame of the Concordance ACDM model. 

PACS numbers: 98.80.-k, 04.50. Kd, 95.36.-l-x 

^ . I. INTRODUCTION 

It is a well-known fact that modifying the law of gravity renders possible explanations for the acceleration mechanism 
of the Universe Q-fl]- However it is far from clear which class of dark energy (DE) theories will finally prevail and 
^ ■ all the viable mechanism must be studied very carefully. Whereas the theories explaining the accelerated expansion 
1^ \ of the Universe in the framework of the general relativity (GR) Q are usually dubbed as DE models, theories in the 
framework of modified gravity are more specifically referred to as modified gravity DE theories. 

In this letter we consider a class of modified gravity theories in which the gravitational action contains a general 
• , function f{R, T), where R and T denote the Ricci scalar and the trace of the energy-momentum tensor, respectively. 
This kind of modified gravity was introduced first by Harko et al. Q where some significant results were obtained. 
In the framework of f{R,T) gravity, some cosmological aspects have been already explored: the reconstruction of 
I ' cosmological solutions, where late-time acceleration is reproduced was studied in Ref. 5] and the energy conditions 
\ analyzed in Ref. j^]. The thermodynamics of Friedmann-Lemaitre- Robertson- Walker (FLRW) spacetimes has been 
studied in Ref. 0, and also the possibility of the occurrence of future singularities (see Ref. Q). 

So far, a serious shortcoming in this kind of theory has been the non-conservation of the energy-momentum tensor. 
\^ . In this paper we circumvent this problem by showing that an algebraic function f{R, T) can always be constructed to 
be consistent with the standard energy-momentum tensor conservation. In the following we shall assume separable 
algebraic functions of the form f{R,T) = fi{R) + f2{T). Within this special choice, the function f2{T) is obtained 
by imposing the conservation of the energy-momentum tensor. 

Once the cosmological background evolution is known, the following step consists of determining the evolution of 
matter perturbations. The growing of structures is manifestly dependent on the gravitational theory under consider- 
ation. This fact can be used to test alternative theories of gravity in order to find out whether those theories are in 
agreement with GR predictions and experimental data 4^ and break the so-called degeneracy problem [T]| that 
some modified gravity theories suffer at the background level. 

This sort of work has been developed in the last years but mainly for the f{R) gravity scenarios 112]. In this realm, 
one of us studied the evolution of scalar cosmological perturbation in the metric formalism within the longitudinal 
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gauge fis'l proving that when scalar cosmological perturbations are studied, f{R) theories, even those mimicking 
the standard cosmological expansion, usually provide a different matter power spectrum from that predicted by the 
ACDM model. 

Still in the framework of f{R) theories of gravity, Gannouji et al. considered linear growth of matter perturbations 
at low redshifts in some DE models and showed the differences with scalar-tensor theories [1^]. For further details 
about cosmological perturbations within f{R) gravity see jisj . 

Therefore the dynamics of cosmological scalar perturbations is a powerful tool to constrain the viability of different 
modified gravity theories, by comparing its density contrast evolution with GR expected features [18, 19]. 

Nonetheless, no attention has been yet paid to study the density contrast evolution in /(i?, T) theories. 

For our purpose in this communication, the dynamics of linear perturbations are performed studying the problem 
of obtaining the exact equation for the evolution of matter density perturbations for fi{R) + f2{T) type gravitational 
Lagrangians. More precisely, we shall assume for simplicity the algebraic function fi{R) as the Einstein-Hilbert term 
R and the trace dependent function f2{T) the one for which the covariant conservation of the energy- momentum is 
accomplished. 

In the so-called quasi-static approximation all the time derivative terms for the Bardeen's potentials are neglected, 
and only time derivatives involving density perturbations are kept 13j,il5i]. Let us finally point out that this approxi- 
mation may remove essential information about the evolution of the first-order perturbed fields [l6l . [22| and therefore 
requires careful study when considered. 

The paper is organized as follows: in Section we briefly review the state-of-the-art of f{R,T) gravity. Section 
mil is devoted to introduce the background cosmological equations for f{R,T) — /i(i?) -I- f2{T) models as well as 
the condition to guarantee standard energy-momentum conservation for such models. Then, Section IIVI addressed 
the calculation of the scalar perturbed equations for f{R, T) = fi{R) + f2{T) models while Section |V] deals with the 
study of the quasi-static approximation for this kind of models. In Section IVll we apply our results to two particular 
models and numerical results are obtained and compared with the ACDM model. Finally in Section [VIII we conclude 
with the main conclusions of this investigation. 



II. f{R, T) GRAVITY THEORIES 

Let us start by writing the general action for /(i?, T) gravities [3], 

S ==SG + S,n = ^ J d^xV^ f{R, T) + J d^x^ , (1) 

where, = SttG, R is the Ricci scalar and T represents the trace of the energy-momentum tensor, T = T^, while 
Cm is the matter Lagrangian. As usual the energy-momentum tensor is defined as, 

2 SSm 

Then, by varying the action with respect to the metric field g'^'^, the field equations are obtained, 

fniR, T)R^, - T)g^, - {g^,U - V^V,) /^(i?, T) = - {k' + MR, T)) T^, - MR, T)Q^, , (3) 

where the subscripts on the function ,f{R, T) mean differentiation with respect to R or T, and the tensor Qf^i, is 
defined as, 

e,. ^ 5"^^ = -2r,. ~ g,^Cm + 2g-l'^^ . (4) 

Note that for a regular f{R, T) function, in absence of any kind of matter, the corresponding f{R) gravity equations 
are recovered, and consequently the corresponding properties and the well-known solutions for f{R) gravity are also 
satisfied by f{R,T) theories in classical vacuum (for a review on f{R) theories, see [l|). Moreover, here we are 
interested to study the behavior of this kind of theories for spatially fiat FLRW spacetimes, which are expressed in 
comoving coordinates by the line element, 

ds^ = 0^(77) (dT]^ -dx^) , (5) 

where 0(77) is the scale factor in conformal time 77. Then, the main issue arises on the content of the Universe is given 
by through the energy- momentum tensor, defined in Since we are interested on flat FLRW cosmologies, the usual 
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content of the Universe (pressureless matter, radiation,... ) can be well described by perfect fluids, which to take the 
form. 



= (P + P)u^^u„ - pg^„ 



(6) 



Here p and p are the energy and pressure densities respectively, and is the four- velocity of the fluid, which satisfies 



1, and in comoving coordinates is given by = (1, 0, 0, 0). Since Cm = P, according to the definition 



suggested in Ref. the tensor (|4]) yields. 



Thus the equations motion become 



(7) 



(8) 



where we have dropped the explicit dependences of / in i? and T . 

It is straightforward to see that the usual continuity equation is not satisfied for the field equations ([8]), and 
consequently the covariant derivative of the energy-momentum tensor is not null in general V^T^" 7^ 0. In order to 
obtain the modified continuity equation, let us take the covariant derivative of the equation ([5]), 



(9) 



Using the identities {Rf_tu — \Rgfiu) — 0, and (V^D — nV^) fR{R,T) = R^p^'^fR, the covariant derivative of the 
energy-momentum tensor needs to satisfy. 



(10) 



Hence, for a perfect fluid with an equation of state p = wp, being w a constant, the 0— component of the covariant 
derivative (fTOl) turns out to become, 



w — 3 



/t - (1 + w)TfTT 



r + 3(l + w) H{K^-fT)-2fTR{AHH + H) 



T = 0, 



(11) 



where let us remind that T = — p — 3p. The last equation differs from the usual continuity equation on the 
non-null right hand side (r.h.s.). T?hus, it may lead to violations of the usual evolution of the different species in the 
Universe. Nevertheless, in the next section we focus our attention on a model that keeps the usual continuity equation 
unchanged. 

III. h{R) 4- /2(T) TYPE THEORIES 

In this section, we choose the algebraic function f{R,T) to be a sum of two independent functions 

/(i?,r) = /i(i?) + /2(r), (12) 

where fi{R) and f2{T), respectively depend on the curvature R and the trace T. The generalized Einstein equations 
from (|5]) yield 



^'RfiRa + S'H'/iflo - y ./lo = -K^a^po + (1 + cl)poa'^f2T„ + y ./20 , 



Aflo + ^fiRo - i^' + 2H2)/iflo + y /lo 



-K a c^po - y/20 



(13) 
(14) 
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where the prime holds for the derivative with respect to rj, Ti = a' /a and the subscript holds for unperturbed 
background quantities: i?o denotes the scalar curvature corresponding to the unperturbed metric, po the unperturbed 
energy density, with /lo = /i(i?o), fiRo = d/i(i?o)/di?o, /20 = /2(7o), /2T0 = d/2(ro)/dro and = Po/po- The 
continuity equation pT|) for Lagrangians given by ((T^) yields 



\7 — 



(15) 



showing explicitly that the energy-momentum tensor is not a priori covariantly conserved in /(i?,T) theories. Thus, 
for these theories, the test particles moving in a gravitational field do not follow geodesic lines. By exploring the 
equation fTS]) for 1/ = component, one gets 



p„ + 3Hpo(l + c?) = 



1 



/2T0 



(1 + C?)Po/2To + ^Ip'ohn + \f'20 



(16) 



Note that whether /2 vanishes (i.e., /(i?) theories) or characterizes a non-running cosmological constant, both f2To 
and /2T0 vanish, and then the continuity equation in these scenarios becomes 



Po + (1 + c2) Po = 



(17) 



In order to Lagrangians such as p2|l consistent with the standard conservation equation, the r.h.s. of (jTS)) has to 
vanish leading to the differential equation 



^l + C^)ro/2ToTo + i(l-c2)/2To=0, 



where 7^ 1/3. The general solution of this differential equation reads 



(18) 



(19) 



where a and /3 are integration constants. In the case of a barotropic equation of state ~ 0, i.e., dust, the model 
(fT9l) becomes 



/2(To)-aTo^/V/?. 



(20) 



This function represents the unique Lagrangian that satisfies the usual continuity equation (|17|) within the class of 
models given by expression (|12p. 



IV. PERTURBATIONS IN f{R,T) THEORIES 



Let us consider the scalar perturbations of a flat FLRW metric in the longitudinal gauge 

ds^ ^ a^[rj) [(1 + 2^)dr]'^ - (1 - 2^)dx^] , 



(21) 



where $ = $(77, x) and 'J = ^'(77,x) are the scalar perturbations. The components of perturbed energy-momentum 
tensor in this gauge are given by 



STl ^Sp^ P06 , ST\ 



-6p S\ 



-cIpoS)5, 5T°, 



-5T\ 



(1 + Cs) Padiv , 



(22) 



where v denotes the potential for the velocity perturbations. Using the model (|12p . the perturbed metric ((2T|) and 
the perturbed energy- momentum tensor ([22|) . the first order perturbed equations reads 



/ifl„<5G;; + (i?ot + V^V, - <5^;n) hn^RjR + '55''" V,V„ - <5£; 



IlRo 



-(1 ~ cDhToK + (1 - 3c?)(l + C^)pof2ToToU''u, 



5p , 



(23) 
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where fiRgRg = d^/i(i?o)/di?o, V"Va and V holds for the covariant derivative with respect to the unperturbed metric 
([5]). In (|23p . we have made use of the relation linking the trace to the enery density, Tq = pg — 3po = (1 — 3c^)poj and 
by the way, ST = (1 — 3c^)Sp. Here the equations of motion at the left hand side of (1^51) presents a set of fourth-order 
differential equations. By following the same assumptions, the equation of the perturbations of the continuity equation 
(fT5|) can be easily obtained, which yields. 



+ Si^d^ { -f2ToST+pof2ToToST + 6pf2To) +dM2nm+TMf2T„T„ST) 



(24) 



from which the conservation equations in /(i?) gravity are obtained when f2{T) is equivalent to the cosmological 
constant, or vanishes. 

For the linearised equation (l23l) . the components (ij), (00), (ii) and (Oi) = (iO), where i,j — 1,2,3, i ^ j in Fourier 
space, read respectively: 



(25) 



{m^ - + W) f[j,„ =a^[- K^po + (1 - 2cl - 3ct) plf2ToTo + ^ (3 - cl) P0/2T0 



(26) 



$" + + m ($' + ■^') + 3-H'$ + {n' + 2n^) * + (sh* -n^ + 3$') f[j^^ 

+ (3$ - VI/) /('^^ = a' [n'dp, + ^ (1 - 3cl) pof2To 



(27) 



(2$ - *) + + + ($ + VET) = -a^ (1 + d) {>,^ - /2To) Po 



(28) 



with 



<5i? = -- 



3*" + 6{n' + n^) $ + 3-H ($' + 3*') -k^{^- 2*) 



(29) 



where it is easy to remark that for /2(To) = the f{R) equations are recovered Moreover, for /i(i?o) = Ro, 

the GR equations are obtained [lO]. Now, by using = 0, ([TSl) and ([24|) . the energy- momentum tensor conservation 
renders to the following first order equations 



S'Hf2ToToPaS ( 1 



(^^-/2To)2 



2/2T0 + P0/2T0T0 



1 



/2T0 



-3H(5 ( -P0/2T0T0 + P0/2T0T0T0 + 2/2T0 



(30) 



and 



—7 \ \hTo5 + 3'Hpof2ToToV 

- /2To V 2 



(31) 



for the temporal and spatial components respectively. 

From the previous expressions is clear that for /2(To) = 0, the usual conservations equations in f{R) theories (GR 
in particular) are recovered (see for instance eqns. (21) and (22) in 13]). Note that expression ([TT]) has been used in 
order to obtain both (f30|) and ([31]) . After further simplifications, the last two expressions become 



S' - k^v - 3*' = 



(32) 



6 



and 



that when combined yield 



5" +% 



if: 



2To 



2(«2-/2To) 



5' + e 



/2T0 

2(«;2-/2To) 



{mv - 5) , 



2To 



2(«:2-/2To) 



1 - 



3/: 



2To 



= 0. 



(33) 



(34) 



Hence, the complete set of equations that describes the general linear perturbations for the kind of models considered 
here, /(i?, T) = fi{R) + f2{T), have been obtained, which provides the enough information about the behavior of the 
perturbations within this class of theories, that can be compare with those results in ACDM model. 



V. EVOLUTION OF SUB-HUBBLE MODES AND THE QUASI-STATIC APPROXIMATION 



We are interested in the possible effects on the density contrast evolution once the perturbations enter the Hubble 
radius in the matter dominated era. In the sub-Hubble limit, i.e., "H <C A:, and after having neglected all the time 
derivative for the Bardeen's potentials $ and ^I*, the equations (l25t and (p6|) can be combined yielding 



1 



2k^ flRpRo 
/iRn 



$ = - 



1 

2fc2 



1 



4fc^ /iHq-Rq 

/iRq 
3fc^ fiRpRo 

1^ flRp 



(«' - /2T0) 



iRp 



In addition, the equation p4[) in the QS approximation yields. 



n 



3/2T0 



2{K^-f2Tp) 



5' + e 



/2T0 



2(«2-/2To) 



5 + = Q. 



Then, by using the previous result (|35l) in the equation ([36|) one gets 



1 - 



3/2T0 



2(«2-/2To) 



/2T0 



(«^-/2To) 



-(«'-/2To) 



/iflo \ 1 + 3 



1 + 44^^^^ 

fc^ flRgRp 



0, 



(35) 



(36) 



(37) 



that can be understood as the quasi-static equation for f{R,T) models of the form ([T^ . By neglecting in (I57|) the 
terms /2('7o), i.e., paying attention only to f{R) theories, one recovers the usual quasi-static approximation for those 
theories (see for instance [2l|, and Q) 



ns 



1 I A fc^ f^RpRp \ 9 , 



1 



1 fe2 /lRo«0 
* — 1 — 2 

1 flRp 



2/iflo 



5 = 



and for GR (/i(i?o) = ^0), the quasi-static equation for S becomes the well-known expression 

S" + US' - A-KGpoa^S = , 



(38) 



(39) 



which is fc-independent. 

Note that the effect of the /2(To) terms in (|37p is twofold: first, the coefficient of 5' gets an extra term that depends 
on the first derivative of /2(To) with respect to Tg that in general will be time dependent. Second, the coefficient 
for 5 is also modified by adding a dependence that is absent the standard quasi-static limit both in GR and in 
f{R) theories and modifying as well the usual coefficient already present for f{R) theories by a factor (k^ — /2To)- 
The fc^-presence may have extraordinary consequences since for /(i?) theories it is usually claimed that in the two 
asymptotic limits (i.e., either GR or f{R) domination), the quasi-static equation is scale independent and only in the 
transient regime, differences associated to the scale may show up. For the class of /(i?, T) theories under study there 
is always a fc^ term that for deep Sub-Hubble modes will be dominant at any time of the cosmological evolution. 

On the other hand, a qualitative analysis taking into account that 
(10~'^eV)"^, implies that equation ([571) may be simplified yielding 



Mp^ « (10i9GeV)-2 and /2T0 



^1/2 

P critical 



S" + -H5' 



l2To- 



fiRo 



' 1 -)- /[fc^ hRpRp 



flRn 



2 _j_ gfc^ fiRpRp 
flRp 



(40) 
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Furthermore, if now one is interested only in extreme sub-Hubble modes, it is clear that pOj) becomes 

= 0, (41) 



2 2 



4 aV) 



that in this limit and after having considered reasonable gravitational Lagrangian, i.e., not divergent, yields 

6" + ^H5' -^k^S = 0. (42) 

The last expression, as well as the intermediate results (j40|) and (|4T|l prove that gravitational Lagrangians depending 

on the trace of the energy-momentum tensor and satisfying the usual conservation equation will exhibit a density 

contrast evolution that is fc-dependent for sub-Hubble modes. This fact is in contradiction with the usual assumptions 

about scale-invariant spectrum of matter perturbations before entering the Hubble horizon that can be seen, for 

instance in the GR density contrast evolution (15^1) . Therefore, models such as (|^D|) . i.e., the only ones guaranteeing 

the standard conservation equation, seem to be theoretically excluded as will be explicitly shown in the next section. 

In addition, note that the equation (jST]) exhibits a singular point at — /2T0 = 0. For the Lagrangian /2(T'o) = 
1/2 

aTg' -|- /3, such singular point is easily identified. From now on, let's assume the following coupling constant, 

a = CiK^po{ritoday)Y^'^ , (43) 

where poiritoday) = T^'^m^o ^^.d Ci is a dimensionless constant, whereas f2^„ and i?o hold for the matter density and 
the Hubble parameter values today (at redshift z = 0) respectively. This parametrization is justified in order to fix 
the correct dimensions for the coupling constant a. On the other hand, by solving the continuity equation p7)) for a 
pressureless fluid, the evolution of the matter density yields, 

PO = Ptodaya^^ = Ptoday(l + , (44) 



where the usual relation 1 + z = a has been used. Then, the expression of the denominator in the equation (|37|) is 
given by, 



2 J- z 
K - /2To = K 



1- 



(l + z)3/2 



(45) 



2/3 

Hence, a singularity occurs at Zs — Ci — 1. Then, the avoidance of such singularity constrains the value of the free 
parameter ci: 

• Cl < 0, the singular point is located at Zg < —1, outside of the allowed range for the redshift, as defined above. 

• Cl > 0, here we can distinguish between two cases: if < ci < 1, then —1 < Zs < 0, and the singularity will 
occur in the future, while if ci > 1, the singularity is located at z^ > 0, so at the present or past cosmological 
evolution. 

In order to avoid any singularity, at least for the range z > 0, we shall assume ci < 1. Note also that in the 
neighbourhood of the singularity, the equation p7p reduces to, 

^" - ^ o, ¥""^1 + I'"^", , S = (46) 

2(k^-/2To) 2(k^-/2To) 

and in consequence the perturbations would behave as a damped oscillator, as is analyzed in the following section 
and shown in Fig. 3. 



VI. NUMERICAL RESULTS 



In order to check the results obtained in the previous section, we study two particular f{R,T) models where 
we have assumed that the /i(i?o) function is represented by the usual term proportional to the Ricci scalar, i.e., 
/i(i?o) = ^0- This choice encapsulates a modification to GR purely originated by the function /2(7o) introduced in 
Section 2 through the expression (fT9|) . 
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Log[1+z] 

Figure 1: 5^. evolution for fA{R,T) model according to the quasi-static evolution given by II40I I and ACDM given by II39I I. The depicted 
modes are k = 50, 100, 200, 500 and 1000 (in Hq units). The plotted redshift ranged from z = 1000 to z = (today). The value of 
was fixed to 0.27 for illustrative purposes. It is seen how whereas the ACDM is fc-independent, the /a(^?, T) model evolutions diverge for 
all the studied modes and leave the linear region at redshifts z 100. For larger fc-modes (deep Sub-Hubblc modes) the divergence 

happens at larger redshift (earlier in the cosmological evolution). 



A. fA{Ro,To)^Ro + aT^^^ 

This model encapsulates a modification to GR purely originated by the function /2(2o) introduced in Section 2 
through the expression ^T9\\ . 

For this function we parametrize the constant a as follows from (|43|) that possesses the appropriate dimensions and 
where ci holds for a non-dimensional constant. For this model, one can solve the background evolution 

^ Ca-i + (l-l]^)aV2 (47) 

numerically where we are using dimensionless time defined as ?7 = Hqi] and f2^„ holds for the usual fractional matter 
density today. According to the last equation the parameter ci must accomplish 

ci = -V^, (48) 

in order to satisfy (|47|) . 

For this model, we compare our result (PH)) with the standard ACDM quasi-static approximation The initial 

conditions are given at redshift z — 1000 where S is assumed to behave as in a matter dominated universe, i.e. 
Skit]) cx a(ry) with no /c-dependence. 

In Fig. 1 we have plotted the evolution of the density contrast for several modes. One can see how the strong 
/c-dependence of equation (|40p renders the evolution of these modes completely incompatible with the density contrast 
evolution provided by the Concordance ACDM model and leads S outside the linear order at redshift z « 100. 



B. fB{Ro,To) = Ro + aTa'^ - 2/3 

Let us now consider the general model found in (|19p . which also satisfies the usual continuity equation in the 
background but where the usual GR term is supplemented with a cosmological constant —2(3. The first FLRW 
equation (fT4|) yields, 

H = 17°„a-i - ciiliy/^ + C2a^ , (49) 
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k=50 H„ 

k=100H„ 

k=200 H„ 




Log[1+z] 

Figure 2: 5^ evolution for fB{R,T) model according to the quasi-static evolution given by II37I I and ACDM given by II39I I. Here we 
have assumed a value ci = —10"^. As previously, the dependence on k leads to a strong growth of the matter perturbations for large 
values of k, whereas the behavior is similar to the ACDM model for the modes k < 200-ffo. 

where we have P5)) . and /3 = iH^C2 in order to provide the correct dimensions to the free constants parameters 
{a, (3}. We have also rewritten the Hubble parameter a.s H = HqH, where Ho is the value today, and have defined 
f2°„ = P"^(^today) ^ gy assuming a{r]today) = 1 (-^ = 0), the equation (|49)) must satisfy, 

i^nl-c,nl + c2 ^ C2 = i-nl,{i-ci) . (50) 

This expression provides a constraint on the dimensionless parameters {ci,C2}, where one remains arbitrary. Thus, 
one can explore the scalar cosmological perturbations at the quasi-static approximation through the equation p7p . As 
for the previous case, the strong dependence on k in the equation p7|) leads to an evolution of the matter perturbations 
incompatible with the observations. In fact, only a very restricted limit for the free parameter ci can avoid such strong 
violations together with an upper limit on k. In Fig. 2, the case for a negative ci = —10^"^ is considered, yielding a 
similar behavior as in Fig. 1. Another illustrative example of the behavior of the equation p7p is shown in Fig. 3, 
where ci = 10""^ is set. For this case, it is straightforward to check that the equation ([57)1 turns out the damped 
oscillator equation for large k — modes, since the k dependent term is positive and dominates over the other terms 
for small redshifts. 



VII. CONCLUSIONS 

In this work we have studied the evolution of matter density perturbations in f{R,T) theories of gravity. We 
have presented the required constraint to be satisfied by these theories in order to guarantee the standard continuity 
equation for the energy-momentum tensor. This constrain restricts severely the form of f{R,T) models able to 
theoretically preserve both BBN abundances and the usual behavior of both radiation and matter as cosmological 
fiuids. 

For models of the form fi{R) + f2{T) we have determined the unique f2{T) cx T^/^ model able to obey the standard 
continuity equation. 

Once the background evolution was imposed such viability condition, we have obtained the quasi-static approxima- 
tion for these theories and shown that for sub-Hubble modes, the density contrast obeys a second order diff'erential 
equation with explicit wavenumber dependence and subsequent strong divergences on the cosmological evolution of 
the perturbations. 

This fact is in contrast with well-known results for f{R) fourth order gravity theories and also Hilbert-Einstein 
action with a cosmological constant. 
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Figure 3: 5^ evolution for T) model according to the quasi-static evolution given by II37I I and ACDM given by II39I I. Here we 

have assumed a positive value for the free parameter ci = 10~^, which leads to an oscillating behavior of the matter perturbations, 
which turns out stronger as k is larger, and whose oscillations are observed for large small redshifts. The model mimics the ACDM 

model only those modes small enough k < 50-ffo ■ 



We have then compared our results with the usual quasi-static approximation in general relativity and shown 
how both density contrast quantities evolve differently. Also the departure from the linear regime happens for these 
theories. These results are in strong contradiction with the usually assumed behavior of the density contrast and 
therefore sets strong limitations to the viability of these theories. 

Alternatively, the study of a positive coupling constant for the modified term T^/"^ leads to a damped harmonic 
oscillator for large fc-modes, as we illustrated in the second model under consideration, in particular in Fig. 3. 
Moreover, as analyzed above, the quasi-static approximation equation contains a singular point that may force the 
matter perturbations to diverge along the cosmological evolution. This assumption provides a way to constraining 
the value of the coupling constant ci, but does not prevent the strong deviation of the sub- Hubble models for this 
kind of models. Hence, we conclude that the kind of models studied in this investigation, f{R,T) = fi{R) + f2{T), 
where the only viable f2{T) function is given by p9|) . leads to a catastrophic behavior for the sub-Hubble modes in 
the study of matter perturbations, preventing this kind of models to be considered as a competitive candidate for 
dark energy. 
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